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1. Introduction and overview 

Nonequilibrium quantum field theory is a research area showing substantially increased 
theoretical activity in recent years. High-energy particle physicists as well as cosmolo- 
gists are working on very similar techniques and sometimes even on the same underly- 
ing nonequilibrium phenomena. This extends to other rapidly evolving topics such as 
nonequilibrium dynamics in laboratory experiments of ultracold quantum gases [|2]. 

Many developments are triggered by high-energy physics related to collision experiments 
of heavy nuclei ("Little Bang") and early universe cosmology ("Big Bang"). Despite sub- 
stantial differences, there are some remarkable parallels. An important example concerns 
the role of nonequilibrium instabilities for the process of thermalization. Though the 
origin of these instabilities can be very different, they always lead to substantial growth 
of occupation numbers in long wavelength modes on time scales much shorter than the 
asymptotic thermal equilibration time. For an anisotropic QCD plasma Weibel instabili- 
ties may operate For the reheating of the early universe after inflation a tachyonic or 
parametric resonance instability can lead to an exponential growth of occupation num- 
bers [Hj. Characteristic far-from-equilibrium phenomena, such as an early prethermal- 
ization of the equation of state, have been quantitatively studied in that context [El E]- 
This is followed by a comparably long quasi-stationary period in a manner reminiscent of 
Kolmogorov wave turbulence [0IE]- The turbulent behavior is described by stationary so- 
lutions of classical statistical field theory. Quantum corrections finally lead to deviations 
and thermalization to Bose-Einstein or Fermi-Dirac distributions [ El EH El E] • 

The presence of strong interactions and/or large occupation numbers from nonequi- 
librium instabilities require nonperturbative approximations. Aspects of systems with 
high occupation numbers can be nonperturbatively described using classical field theory 
methods. However, classical Rayleigh- Jeans divergences and the lack of genuine quan- 
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turn effects — such as the approach to quantum thermal equilibrium characterized by 
Bose-Einstein or Fermi-Dirac statistics — limit their use. Gauge-string duality offers a 
novel possibility for studying dynamical properties of certain strongly interacting gauge 
theories, which has led to interesting results [IT5]. 

Standard nonperturbative approximations such as based on 1/N expansions of the 
one-particle irreducible (1PI) effective action can be secular in time and do not provide 
a valid description, similar to perturbation theory. Suitable resummation techniques can 
be efficiently formulated using the two-particle irreducible (2PI) effective action [ 14J. 
This has led to successful descriptions of far-from-equilibrium dynamics and subsequent 
thermalization in scalar and fermionic quantum field theories in various dimensions [ 
ITT] ITB1 Uni El El Ell U21 El ■ They have also been used to compute transport coefficients [ 
ITH] and to determine the range of validity of transport or semi-classical approaches [ 

nulla eh 

A nonperturbative description of nonequilibrium quantum fields can be based on the 2PI 
1/N expansion beyond leading order (LO), which has been worked out in detail [rHHl2"2lE)l 
EJEHEIE]- S° far, the results from the 2PI 1/N expansion to NLO are the only ones that 
bridge nonequilibrium instabilities at early times and thermal equilibrium at late times. 
Recently the 2PI 1/N expansion has been pushed forward towards NNLO for a quantum 
anharmonic oscillator [I2H]- A series of precision tests has been performed by now [1231213 
ITU] , including the computation of critical exponents near second-order phase transitions [ 
I2H1- The results exhibit some remarkable properties of the 2PI-l/iV expansion, curing for 
instance the spurious small- N divergence of the anomalous dimension seen in the standard 
1PI 1/N to NLO. Recently, also the question of the formation of topological defects has 
been investigated in the large- N expansion for N = 1, 2 in one spatial dimension at zero 
temperature. Since there are no defects in the corresponding quantum theory at zero 
temperature, the authors drop quantum corrections and introduce a cooling procedure. 
Defects should then appear as an artefact of the classical approximation, however, the 
authors see no sign of them. The size of the observed discrepancies are in accordance with 
earlier comparisons within classical statistical field theory [[2Zj- An important source of 
quantitative differences arises from the fact that a large- N expansion to NLO in the limit 
N — > 1 trivially misses 1/3 of the topologically equivalent diagrams starting at NNLO [ 

1221- 

A nonperturbative large- N expansion of the 2PI effective action beyond LO, feasible for 
the (abelian) case of large numbers of flavors, seems not within reach for SU(N) gauge 
theories relevant for QCD. In this context it remains an important open question to 
determine the range of validity of 2PI loop expansions, in particular, in view of a number 
of remarkable results from related Dyson- Schwinger equations in thermal equilibrium [ 
I2H]. It has been an important step that renormalization of 2PI resummed approximations 
for gauge theories has been worked out [ I2H1 EOj , extending earlier work for scalar [ EJJ 
and fermionic theories [I32j. 

Complementary to approaches based on analytic approximations, one of the outstanding 
problems in nonequilibrium quantum field theory are first-principles simulations on a 
Minkowski space-time lattice. In the following we will describe recent efforts along these 
lines, which are still in its infancies. 




Figure 1. The real part of (J (f(0,x)f(t,x)) as a funtion of real time t in units of the 
lattice spacing. Shown are snapshots at different Langevin times The LHS corre- 
sponds to free field theory, with null starting configuration at $ = 0. The RHS is for the 
interacting theory. Here the nonequilibrium classical statistical result is taken as starting 
configuration, and the Langevin updating for $ > incorporates quantum corrections [ 

EZj. 



2. Lattice simulations of real-time quantum fields 

Lattice gauge theory calculations are typically based on a Euclidean formulation, where 
the time variable is analytically continued to imaginary values. By this the quantum the- 
ory is mapped onto a statistical mechanics problem, which can be simulated by importance 
sampling techniques. Recovering real-time properties from the Euclidean formulation is a 
formidable problem. Direct simulations in Minkowski space-time would be a breakthrough 
in our efforts to resolve pressing questions, such as early thermalization or the origin of 
seemingly perfect fluidity in a QCD plasma at RHIC. For real times standard importance 
sampling is not possible because of a non-positive definite probability measure. Efforts 
to circumvent this problem include mimicking the real-time dynamics by computer-time 
evolution in Euclidean lattice simulations [EHUS]- A problem in this case is to calibrate 
the computer time independently of the algorithm. In principle, direct simulations in 
Minkowski space-time may be obtained using stochastic quantization techniques, which 
are not based on a probability interpretation [ |33J EH] . 

In Ref. [ EZj this has been recently used to explore nonequilibrium dynamics of an 
interacting scalar quantum field theory. As an example, Fig. ^ shows the time evolution 
for the connected part of the unequal-time correlator Re( j <p(0, x)ip(t, x)), which measures 
the correlation of the field <p at time t with the initial field. Shown are snapshots of the 
correlator at different Langevin-times. The left figure shows the result for the free field 
theory with null start configuration. The Langevin updating reproduces the correct free- 
field result. The right figure gives an example for the interacting theory, which exhibits 
a finite characteristic damping time. One observes good convergence properties of the 
quantum simulations, which is a remarkable result. For given initial field configurations 
at time t = 0, very different starting configurations for the (3 + l)-dimensional space-time 
lattice converge to the same nonequilibrium dynamics for all t > [13*7]. 
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2.1. Real-time stochastic quantization 

In real-time stochastic quantization the quantum ensemble is constructed by a stochas- 
tic process in an additional "Langevin-time" using the reformulation for the Minkowskian 
path integral [|3BJIIIII| : The quantum fields are defined on a d-dimensional physical space- 
time lattice, while the updating procedure employs a Langevin equation with a complex 
driving force in an additional, unphysical "time" direction. Though more or less formal 
proofs of equivalence of the stochastic approach and the path integral formulation have 
been given for Minkowski space-time, not much is known about the general convergence 
properties and its reliability beyond free-field theory or simple toy models [ EH] • More 
advanced applications concern simulations in Euclidean space-time with non-real actions [ 

mum. 

Below we follow Ref. [H2| and discuss real-time stochastic quantization for a scalar the- 
ory and SU(N) pure gauge theory relevant for QCD. For a scalar theory with Minkowski 
action S[ip] the Langevin updating equation reads 

(p\x) = (p(x)+ie^-j^r + y/eri(x), (1) 
d(p(x) 

with Gaussian noise 

Ms)),, = , (n(x) n{y)) v =26(x-y). (2) 

The sum over all Langevin steps, *d = ^ e, corresponds to Langevin-time. 

For SU(N) gauge theory on a (N s a s ) 3 x N t a t lattice the real-time classical action reads 

S[U] = -AEEl^I^^ + Tr^j-l} 

E E { 2ttt (Tr Ux > ij + Tr u ^ ~ 1 } ' 



(3) 



1,3 



with spacial indices i, j = 1, 2, 3. It is described in terms of the gauge invariant plaquette 
variable 

U x ^ u = U X) pU x +fi )V U x+i)>IJI U x>l , . (4) 

Here U x ^ is the parallel transporter associated with the link from the neighbouring lattice 
point x + p, to the point x in the direction of the lattice axis \x = 0, 1, 2, 3. The couplings 
are 

R - 2 ^ Trl o - 2Trl 

P0 = n ' ^* = — 9 ' \° 

% 9sl 

with the anisotropy parameter 7 = a s /a t . The Langevin updating equation for U x>tl then 
reads 

U' X41 = exp {i\ a (eiD x ^ a S[U] + \fer] XIM ) } U X)IX , (6) 
with 

1 3 

iD XIM S[TJ] = ~—J2^ T H X aUx,,C x ^ u -C x ^M-lXa). (7) 

i/=0 
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For a compact notation we have defined = (3 S and (3oi = Ao = —A) and 
r7 tt - zt -1 zr^ 1 + zt -1 zt -1 zr . 

C x ,iiv — U X)U U X +£,^U X+ ^^ U + U x _p v U x —u^U x j r jx—o,v ■ (8) 

It is shown in Ref. [ 132] (see also Ref. [ESI) that stationary solutions of equations (JTJ) 
or © always fulfill the infinite set of (symmetrized) Dyson- Schwinger identities of the 
respective quantum field theory. For instance, for the scalar theory one finds from the 
fixed points of the Langevin equation (JTJ the infinite hierarchy of identities 



8S[p] \ 
6<p(x)/ 



= 0, (9) 



2i ((ip(x)) v 8{y -z) + (<p(y)) v 5(x - z) + (<p(z)) v 5(x - y)) , (11) 

and correspondingly for the higher n-point functions in Minkowski space-time. If the 
Langevin updating converges, which is the case for the implementation of Ref. [E2J, then 
it solves the correct infinite set of Dyson-Schwinger identities. In addition, Euclidean 
theories with real actions can be shown to have a unique solution based on positivity 
arguments. A similar argument fails for real-time stochastic quantization. In general, here 
the correct fixed point cannot be chosen a priori without implementing further constraints. 

2.2. Precision tests for scalar and nonabelian gauge theories 

In Ref. [ H2J tests for real-time stochastic quantization are discussed. Similar to what 
has been observed in Ref. [ |37] for scalar fields, also for pure gauge theory one finds that 
previously reported unstable dynamics represents no problem in practice: A combination 
of sufficiently small Langevin step size and the use of a "tilted" real-time contour leads 
to converging results in general. The employed procedure respects gauge invariance and 
appears to be well under control. This is exemplified for SU(2) gauge theory in 3 + 1 
dimensions and for a scalar theory in zero spatial dimension, i.e. a quantum anharmonic 
oscillator. For the latter the Schrodinger equation can be solved as well numerically by 
diagonalization of the Hamiltonian, for comparison. 

For the scalar theory one finds that stochastic quantization accurately describes the 
time evolution for lattices with sufficiently small real-time extent. This concerns nonequi- 
librium and equilibrium simulations at weak as well as strong couplings. As an example, 
Fig. El shows results for the time evolution for 4! scalar theory in zero spatial dimen- 
sion. Shown is the two-point correlation function (</?(0) </?(£) as a function of real time 
t. The simulations are done in thermal equilibrium, in which case the real-time contour 
extends along the imaginary axis to the inverse temperature (3. For the left figure the 
real-time extent of the contour is At = 0.5. Here the upper branch of the contour has a tilt 
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Figure 2. The unequal-time two-point correlation function for the quantum anharmonic 
oscillator as a function of real time t. The (short-time) results obtained from stochastic 
quantization agree to very good accuracy to those obtained from directly solving the 
Schrodinger equation. For the left figure a strong coupling A = 96 and (3 — 1 is employed, 
whereas A = 6 and (3 = 8 for the right. [H2] 



of 0.001/5, so it is almost horizontal and, therefore, realizes to high accuracy a real-time 
contour. For comparison to the stochastic quantization result the time evolution from 
the Schrodinger equation in Minkowski time as well as using the employed complex times 
from the corresponding matrix algebra is given. They all agree to very good accuracy. 

However, when the real-time extent of the lattice is enlarged the stochastic updating 
does not converge to the correct solution. Fig. El shows the two-point correlation function 
as a function of real time t for A = 24 and (3=1. For a short real-time lattice corre- 
sponding to tfinai = 0.8 one observes from the left figure excellent agreement of stochastic 
quantization and Schrodinger equation results. However, doubling the extent of the con- 
tour leads to a qualitatively different behavior as given in the right plot of Fig. This 
difference persists also on finer real-time grids. The non-vanishing imaginary part of the 
equal-time correlator and the loss of time-translation invariance reflects a non-unitary 
time evolution. 

Doing the equivalent investigation for SU(2) pure gauge theory for thermal equilibrium, 
one finds a similar behavior with an important difference. The left of Fig. 0] shows the 
Langevin-time evolution of the spatial plaquette average. This quantity is computed for 
different complex contours, with the approximation that g = g s . The solid line shows 
the result for vanishing real-time extent, i.e. for a Euclidean contour corresponding to 
an inverse temperature (3 = 4. The different dashed curves correspond to results for 
complex contours on isosceles triangles each having a different tilt a with respect to the 
real-time axis. Here a = would correspond to an infinite extent along the real axis. 
One observes that with increased real-time extent or smaller tilt a the correct thermal 
solution is approached less accurately. In particular, it is only approached at intermediate 
Langevin-times, irrespective of the details of the non-vanishing real-part of the contour. 
This aspect differs from the scalar case where short real-time extents lead to stable thermal 
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Figure 3. Real and imaginary part of the two-point correlation function for the quan- 
tum anharmonic oscillator as a function of real time t. Left: for short real-time extent 
of the lattice the stochastic quantization results agree well to those obtained from the 
Schrodinger equation. Right: Equal-time correlator as a function of time (index of lattice 
site along the time-contour). Compared are two simulations where the real-time extent of 
the lattice differs by a factor of two. The larger lattice leads to a qualitatively different, 
non- unitary behavior. [ W2\ 



solutions. 

As mentioned above, all converging solutions of the stochastic dynamics method fulfill 
the same infinite set of (symmetrized) Dyson- Schwinger identities of the quantum field 
theory. This is remarkable in view of the different "physical" and "unphysical" solutions 
that are observed. In Fig. |S] this is visualized for the example of the Dyson-Schwinger 
equation for a spatial plaquette variable in SU(2) gauge theory. Plotted are separately the 
LHS and the RHS of the Dyson-Schwinger equation displayed graphically on the right of 
Fig. The left plot displays the respective LHS and RHS as a function of Langevin-time 
for a typical contour with tilt < a < tt/2. The flow with Langevin-time quickly leads 
to a rather accurate agreement of both sides such that the Dyson-Schwinger equation is 
fulfilled. However, after some Langevin-time they start deviating again, finally leading to 
another stationary value where the LHS and RHS agree to reasonable accuracy. 

The evolution with Langevin-time is governed by fixed points and, for the typical case 
considered in Fig. the correct thermal fixed point is approached at first. However, it is 
not stable and the Langevin flow exhibits a crossover to another (stable) fixed point. 

In conclusion, short real-time physics in thermal equilibrium can be reproduced rea- 
sonably well if the length of the real time contour is small on the scale of the inverse 
temperature j3. For longer contours the boundary conditions in physical time do not 
seem to constrain enough the Langevin flow and the 'life-time' of the thermal (physical) 
fixed point decreases (nonabelian gauge theory), or the fixed point becomes fully unsta- 
ble (scalar theory). A second, apparently stable fixed point develops for large Langevin 
times. This latter represents an unphysical, non-unitary regime, to be recognized by non- 
translational invariant expectation functions and violation of the unsymmetrized Dyson- 
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Figure 4. The spatial plaquette average for SU (2) gauge theory in thermal equilibrium as 
a function of Langevin time. Shown are results for different complex contours. Here a = 
corresponds to a contour with an infinite extent along the real-time axis, while a = tt/2 
denotes the Euclidean contour. The longer the real-time component of the contour the 
less accurate the thermal solution is approached. Here ^crossover denotes the Langevin-time 
at which the crossover to the non-unitary fixed point occurs. [ 14*2"] 



Schwinger identities (while the symmetrized ones are still satisfied, indicating conver- 
gence) [E2j- Further work has to concentrate on further means for controlling and opti- 
mizing the method. The method allows for quite some flexibility as to which quantities 
are chosen to define the stochastic process, or introducing a stochastic re-weighting. 
We are indepted to D. Sexty and I-O. Stamatescu for collaboration on this topic. 
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